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THEORETICAL AND EXPERIMENTAL INVESTIGATION OF HEAT 
TRANSFER BY LAMINAR NATURAL CONVECTION 
BETWEEN PARALLEL PLATES ‘ 

By A. F. Lietzke 


SUMMARY 

Results are 'presented of a theoretical and experimental 
investigation of heat transfer involving laminar natural con- 
vection of fluids enclosed het'ween parallel “walls oriented in the 
direction of the body force, “where one “wall is heated uniformly, 
and the other is cooled uniformly. For the experimental “work, 
parallel walls were simulated by using an annulus “with an 
inner-to-outer diameter ratio near 1. 

The results of the theoretical investigation are presented in the 
form of equMtions for the velocity and temperature profiles and 
the ratio of actual temperature drop across the fluid to the 
temperature drop for pure conduction. No experimental 
measurements were made of the velocity and temperature profiles, 
but the experimental results are compared with theory on the 
basis of the ratio of the actual temperature drop to the temperature 
drop for pure conduction. Good agreement was obtained, 
between theory and experiment for axial temperature gradients of 
40 ° F per foot or larger. 

INTRODUCTION 

Increased application of heat transfer to and from fluids in 
channels has recently required further knowledge as to the 
heat-transfer coefficients and temperature profiles occurring 
with natural convection. Turbine-blade and nuclear-reactor 
cooling are two of the fields concerned with this problem. 
Work on free-convection heat transfer over a vertical plate 
gave good agreement between theory and experiment. Few 
results have been obtained for the similar case of flow in 
channels. Reference 1 obtains a theoretical solution for 
free-convection heat transfer for fluids enclosed in channels. 
The reference uses a postulated velocity distribution to obtain 
a solution. References 2 and 3 extend this analytical work to 
give an exact solution of the equations in more general form 
for constant wall temperature and constant heat flux, 
respectively. Reference 3 includes the effect of forced as 
well as natural convection. These three references treat the 
case of infinite channels with the channel axis oriented in the 
direction of the body force and are subject to the same 
assumptions ; namely, two-dimensional laminar flow, uniform 
axial temperature gradient, and constant fluid properties, 
except that the density is allowed to vary in the bouyancy 
term. 

The purpose of this report is to compare the results of 
theory and experiment as a check on the assmnptions in- 


volved in the analytical work. For mathematical simplicity, 
the flow between two infinite parallel plates was considered 
with one plate heated uniformly and ’the other cooled uni- 
formly. An exact solution was obtained with the assump- 
tions cited previously for references 1 to 3. In order to 
simulate the case of infinite parallel plates experimentally, 
an annulus formed by two concentric tubes was used. For 
all practical purposes, the walls of these tubes can be con- 
sidered parallel if the ratio of the radii of the tubes is near 1 . 
This radius ratio limits the spacing between the confining 
walls (and, consequently, the range of Grashof numbers 
obtainable) unless very large diameters are used. 

A comparison between theory and experiment is given in 
this report on a heat-transfer basis alone, since no temper- 
ature or velocity profile measurements were made. The 
work was done at the NACA Lewis laboratory. 


ANALYSIS 

Steady-state heat transfer through a fluid enclosed by two 
infinite parallel plates oriented in the direction of the body 
force is considered. One plate is heated uniformly and the 
other is cooled uniformly. The flow is laminar and parallel 
to the body force. The fluid properties are assumed con- 
stant, except that the density is allowed to vary in the 
buoyancy term. Viscous dissipation and work against the 
force field are neglected. It is further assumed that the 
axial temperature gradient is a constant throughout the 
system for any particular set of conditions. 

When the aforementioned conditions are applied, the 
energy equation reduces to 


k dH _ dt 
pCp ^ bz 


(Bl) 


(Symbols are defined in appendix A, and a detailed dis- 
cussion of the analysis is given in appendix B.) 

Similarly, the Navier-Stokes’ equations reduce to one 
equation in the s-direction; , 




(B7) 


The reference temperature tf is taken at the center of the 


channel where the net viscous force is zero 


fb^w 


0 


■ Supersedes NACA TN 3328, "Theoretical and Experimental Investigation of Heat Transfer by Laminar Natural Convection Between Parallel Plates,” by A. F. Lietse, 1951. 
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For a closed system of unit width with zero net mass 
through-flow and constant density, 


j: 


wdx=0 


where x is taken to be zero at the center of the channel. 

From these three equations and the appropriate boundary 
conditions 


=_2l' 

\i>X/x^r k 


( 


^xj x= —T 


V^X=-T = 0 


the equations for the velocity and temperature profiles in 
dimensionless form are, respectively, 


W* = 4:V 


r=- 


2v 


( 


/ . VX vx 

1 sm — cosh — 

. 1 VX 

sinh — 

vx \ 
cos — \ 

1 

r 

r 

r 1 

1 sin V 

cosh V 

sinh « 

COS V 1 

\sinh v' cos 2 ) 

sin V 

cosh V/ 


(B17) 


vx 


vx 


sinh V cos u-fsinh — cos 
r r 


cosh V . sin 
cos V ~sinh v 


VX . vx 

cosh V sin w-j-cosh — sin 


sinh V cos v 
sin V "‘cosh v 


) 


(B19) 


where, by definition, 


v=^- 


TrGr 


64 


Py. 

k 


Gr- 


■ 


The temperature drop across the fluid from wall to wall 
is then found from equation (B19) when x—r: 

sin^ i)+sinh^ v 


kv \si 


sinh v cosh i;+sin v cos 


i) 


(B20) 


The temperature difference across the fluid for pure con- 
duction is 


{er)c=~ 


k 


(B21) 


The ratio of the actual temperature drop to that for pure 
conduction is, therefore, 

sin^z)+sinh®« 




sinh V cosh sin v cos v. 


) 


(B22) 


These equations apply equally well to force fields other 
than gravitational, if the constant g is replaced by the con- 
stant of the force field of interest. 

EXPERIMENTAL APPARATUS 

A schematic diagram of the experimental equipment is 
shown in figure 1. In order to simulate infinite parallel 


plates, an annulus formed by two concentric stainless-steel 
tubes was used to contain the fluid. The resulting annulus 
was 10>2 inches long with an outside diameter of 1% inches 
and an inside diameter of 1 inch. Thus, the spacing s be- 
tween the walls was inch. The tube walls were Xe inch 
thick. The outer tube was heated with a Nichrom e-wire 
element spiral-wound around the tube, with a wire spacing 
of %2 inch, which probably gives fairly uniform heating of the 
hot wall. The inner tube was cooled by forced air. The 
cooling-air-passage length-to-diameter ratio was 84. A wire 
in the form of a helical spring was inserted in the cooling 
passage to increase the heat-transfer coefficient. This 
system should provide an essentially constant heat-transfer 
coefficient, which is necessary to provide uniform cooling of 
the cold wall with a linear variation of wall temperature. 
The cooling-air-flow rate was measured with a calibrated 
rotameter. 



Figure 1. — Schematic diagram of experimental apparatus for measuring 
heat transfer by free convection in an annulus used to simulate 
parallel plates. 
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Figure 2. — Dimensionless velocity profiles for various values of product 
of Prandtl number and modified Grashof number. 



Figure 5. — Temperature measurements of heated and cooled walls of 
test section. Heat flux, 1300 Btu per hour per square foot; Prandtl 
times Grashof number, 3480; temperature-differance ratio, 0.389. 


measured with iron-constantan thermocouples and a self- 
balancing potentiometer. The thermocouples used to meas- 
ure wall temperatures were set into slots milled axially into 
the tube walls. The thermocouple leads were enclosed in 
Inconel tubing laid in the slots. The thermocouple junc- 
tions were covered with insulating cement, and the slots were 
then filled with silver solder and finished off flush with the 
tube surface. 



Figure 3. — Dimensionless temperature profiles for various values of 
product of Prandtl number and modified Grashof number. 



Figure 4. — Theoretical variation of ratio of actual temperature differ- 
ence across fluid to temperature difference for pure conduction with 
product of Prandtl number and modified Grashof number. 


RESULTS AND DISCUSSION 
ANALYTICAL 

The dimensionless velocity and temperature, which can 
be calculated from equations (B17) and (B19), respectively, 
are shown graphically in figures 2 and 3 for a few values of the 
product of Prandtl and modified Grashof numbers. The 
equations show that the dimensionless velocity and tempera- 
ture depend only on this product and the position in the 
channel. The Grashof number resulting from the analysis 
differs from the conventional Grashof number insofar as the 
temperature difference usually appearing is replaced by the 
temperature gradient and the characteristic length appears 
as the fourth power instead of the cube. The product of 
Prandtl number and Grashof number is sometimes referred 
to as the Rayleigh number. 

The velocity profiles of figure 2 show the point of maxi- 
mum velocity moving closer to the wall as the Rayleigh num- 
ber is increased. It is felt that this fact, in combination with 
the stabilizing influence of the wall, would tend to inhibit 
the occurrence of turbulent flow. 

The contribution of convective flow in the heat-transfer 
process can best be seen from the ratio of the actual tempera- 
ture drop across the channel to the temperature drop for 
pure conduction. This ratio, which is the same as the ratio 
of the molecular to the apparent conductivity, can be ob- 
tained from equation (B22) and is plotted in figure 4. For 
large values of PrGr, the actual temperature drop across the 
fluid is less than K the temperature drop for pure conduction. 

EXPERIMENTAL 


The electrical power was supplied by a transformer with a 
variable secondary voltage and rated at 1 kva at maximum 
voltage. The electrical power input was measured with a 
calibrated ammeter and voPmeter. All temperatures were 

348261 — 56 2 


A typical plot of measured wall temperatures for large 
temperature gradients is shown in figure 5. With uniform 
heating of the hot wall and constant heat-transfer coeflBcient 
on the cold wall, the temperature gradients on the hot and 

- '-...T- -j -j; •■j'"'--- . 
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cold walls are constant and equal, except for the extreme ends 
of the test section. Data were not obtained to verify the 
assumption of uniform temperature gradient throughout the 
fluid. 

For small temperature gradients the wall temperatures 
were linear, but the gradient on the heated wall was larger 
than that on the cooled wall. The assumption of imiform 
temperature gradient throughout the system is, therefore, 
invalid for small temperature gradients. With water as the 
test fluid and with e-inch spacing between the plates, the 
minimum value for which the assumption of uniform tem- 
perature gradient is valid is approximately 40° F per foot. 
Additional data are required to determine this limit more 
accurately. The data presented in this report are for 
temperature gradients of 40° F per foot or greater. 

Because of the end effects on the flow and because of the 
heat losses from the ends of the test section, it was impos- 
sible to obtain a heat balance between the electrical heat 
input and the enthalpy rise of the coolant. Measurements 
of the radial temperature drop through the insulation sur- 
rounding the heater indicated a negligible radial heat loss. 
In that section of the tube where the wall temperature varies 
linearly, there is no net axial conduction in the wall. Hence, 
the uniform heat flux (as calculated by the electrical heat 
input) that enters the wall is transmitted directly to the 
fluid. The heat loss is limited to the ends of the test sec- 
tion and, therefore, does not affect the results of the test. 

The calculated temperature drop through the walls con- 
taining the water was negligible compared with the tempera- 
ture drop across the water, and, therefore, the measured wall 
temperatures were taken to be the temperature of the sur- 
face adjoining the water. 

The temperature drop across the fluid for pure conduction 
Ate calculated from the equation for an annulus given 
in reference 4 : 

.4 loge(r„/ri) 

where g'' was calculated from measurements of the electrical 
heat input. 

The physical properties used in evaluating Prandtl num- 
ber, Grashof number, and the temperature drop for pure con- 
duction were obtained from reference 5 for saturated liquid 
water, except for the coefficient of thermal expansion fl. 
Values of the coefficient of thermal expansion of water were 
taken from reference 6 and are plotted in figure 6. These 
values are mean values for 1° temperature change. Inas- 
much as /3 is a function of temperature, different mean 
values would be obtained for larger temperature increments. 
The properties were all evaluated at the bulk temperature 
of the water for each test. For flat plates, the bulk tem- 
perature is the center-line temperature halfway between the 
ends. 

Figure 7 shows the ratio of the actual temperature drop 
across the fluid to the temperature drop for pure conduction 
plotted against the product of Prandtl and Grashof numbers. 


The data points on the figure are the experimental results 
taken with water as the test fluid, and the curve is the result 
of analysis taken from figure 4. 

COMPAKISON BETWEEN ANALYSIS AND EXPEBIMENT 

Inasmuch as no temperature measurements were made 
within the fluid, a comparison between theory and experi- 
ment can be made only for the over-all heat- transfer results. 
Except for the random scatter of the heat-transfer data shown 
in figure 7, good agreement is obtained between theory and 
experiment. Inasmuch as the same assumptions were used 
in references 2 and 3, the agreement shown here lends sup- 
port to the more general analytical treatment given in these 
references. 



Figure 6. — Variation of coefficient of thermal e.'cpansion of water with 
temperature for 1° change in temperature. 



Figure 7. — ^Variation of ratio of actual temperature difference across 
fluid to temperature difference for pure conduction with product of 
Prandtl number and modified Grashof number. 
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The effect of the curvature of the walls on the experi- 
mental results is unknown; but this effect is believed to be 
mhiimized, because the data are in the form of the tempera- 
ture ratio with this effect included in both the numerator 
and the denominator. A theoretical investigation for an 
annulus is necessary to permit a comparison between theory 
and experiment with larger diameter ratios for the annulus. 

In order to determine the upper limit for which the anal- 
ysis is valid, it is necessary to obtain data at higher values 
of the product of Prandtl number and modified Grashof 
number. 
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CONCLUSIONS 

For the range of conditions investigated, the simplifying 
assumptions used in the analysis for free convection in 
channels oriented in the direction of the body force are 
reasonable and lead to accurate quantitative results. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 8, 1954 


APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 
Ai, A 2 . . . constants 


Bi, B 2 . . . constants 
C constant of integration 

Cp specific heat at constant pressure, Btu/(lb)(°F) 

D operator, d/dx 

Fj body force in g-direction, lb/(hr2)(sq ft) 


Gr 

g 

i 

k 

Pr 

P 

<L'o 

r 

Vi 

To 

s 

t 


modified Grashof number, , > dimensionless 

M 

acceleration due to gravity, 4.17X10® ft/hr^ 

V=T 


thermal conductivity, Btu/(hr)(sq ft)(°F/ft) 

Prandtl number, Cp^^|k, dimensionless 

pressure, Ib/sq ft 

heat flux at wall, Btu/(hr)(sq ft) 

heat flux at outer wall of annulus, Btu/(hr)(sq ft) 

1/2 plate spacing, s/2 

inside radius of annulus, ft 

outside radius of annulus, ft 

plate spacing, ft 

temperature, °F 


h 

tr 

t-r 

U 

V 

V) 


w* 

X 

z 

8 

0 

Or 

(0r)c 

e* 

P 

Vt 


reference temperature, °F 
temperature at x=r, °F 
temperature at x— —r, °F 
axl^/2, dimensionless 


, ^ * iPrGr 

velocity, ft/hr 


at 

WpCpS ^ 

dimensionless velocity, ^77 

transverse coordinate, ft 
longitudinal coordinate, ft 
coefficient of thermal expansion, °F~* 


°F 

tr-t-r, °F 


temperature difference across fluid for pure con- 
duction 


dimensionless temperature, 


(t — t^r)k 
q"s 


viscosity, lb/(hr)(ft) 
density, Ib/cu ft 

ratio of actual temperature difference across fluid 
to temperature difference for pure conduction 


APPENDIX B 


DETAILS OF ANALYSIS 


For the conditions of the problem stated in the ANALY- 
SIS, the general energy equation reduces to 


pCp ^ bz 


(Bl) 


Similarly, the Navier-Stokes’ equations for incompressible 
flow reduce to 

bp 


bx 


=0 


bp 

bz" 




b^w 

bx^ 


The body force under the action of gravity is 

Fz = —gp 


(B2) 


(B3) 


Since the pressure is independent of x, then ^ is inde- 


pendent of X and can be evaluated at any value of x. It is 

convenient, however, to evaluate ^ where ^-^=0. Then, 
’ ’ bz bx^ 

from equations (B2) and (B3), 

bp_ 
bg" 


-pf 


(B4) 


where the subscript j denotes a reference condition where 
^^=0. Combining equations (B2), (B3), and (B4) gives 


gp 


0-7)- 


b^w 


(B5) 


From the definition of the coefficient of thermal expansion /3, 


^=1+8 (f-t,) 
P 


(B6) 
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Combining (B5) and (B6) yields 


S=f 


(B7) 


The reference temperature f/ is a function of z and varies 
linearly with z as does t, according to the assumption of 
constant temperature gradient. The velocity w is therefore 
only a function of x, and its derivative is written as a total 
derivative. 

Differentiating equation (B7) with respect to x gives 


d?w gpfi bt 

di? p. hx 

Differentiating with respect to x again gives 

d*w_ gpfi hH 

dx^ p 

Combining equations (Bl) and (B9) yields 
d*w_ gp^fiCj, bt 


dx'^ 


Leta^=^^|-^- Then 
pk bz 


d*w 

dx* 


pk bz 


Ta^=0 


w 


(B8) 


(B9) 


(BIO) 


The solution of this homogeneous equation can be obtained 
by using the operator D : 

(D^+a^)w=0 

{D-\-i^[la){D~i■^|la){D-^^i a)(D~^ia)w^O 

\ ^ 'I 1 

but and Therefore, 


{n. 


V2 




w=0 (Bll) 


w=0 




dw , ?■— 1 _ 

w ^2 

log w+— ax=C 

V2 


i-1 


V2 = 


W 


i— 1 _1°£ 

d'i =Aie V2gV2 


(B12) 


Similarly, from equation (Bll), 


j^iax ax 

w=A2e ■^‘^e 


iax ax 

w—A^e 

iax ^ 
w=Aie'^ 


(B13) 

(B14) 

(B15) 


The general solution is equal to the sum of equations (B12) 
to (Bl5), or 

w=^e-^^(Aie^+A^e-’‘)+e‘^(A2e-’‘+A,e^) 

, ax 

where u——p-- 

V2 

Using the Euler formulas 


^-iu. 


= cos u—% sm u 


gii/_cos u+f sin u 
and the conversion formulas 

e“=cosh w + sinh u 
g-“=cosh sinh u 

and letting 

Bi=^ A\-\- Ai-\- Ai-\- Ai 
B2 —AiA Ai — A3 A2 
B3=i{A2 — AsAAi — Ai) 

Bi=i{Ai — Ai — A2~\-A3) 

the velocity equation becomes 

w=Bi cos u cosh U + B 2 cos u sinh u + B^ sin u cosh U + B 4 

sin u sinh u (BI6) 

The constants of equation (B16) can be evaluated from 
the boundary conditions. There are four boundary condi- 
tions required to determine the four constants. The velocity 
at the bounding walls must be zero, while the boundar}’’ con- 
ditions on the temperature are given by the temperature 
gradients. It can be shown from equation (Bl) that, for 
the case of interest here, the temperature gradients at the 
two walls must be equal. Integrating equation (Bl) gives 


PCj,Wbx)i:=r \C)3:A=-rJ C>2j_ 


w dx 


For a system closed at both ends, there is no net mass 
through-flow; therefore, with constant density. 


J wdx=b 
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Consequently, the transverse temperature gradients at the 
two walls must be equal. The boundary conditions are 
expressed mathematically as 


\Qa:A=r k 

(^\ =_2! 
, k 


Wx_r=0 

Wo:=_, = 0 

When these boundary conditions are imposed on equations 
(B8) and (Bl6), the constants of equation (Bl6) are found 
to be 

Si=0 

Aq^'v sin v cosh v 




B, 


(sin V cos v+cosh v sinh v) 




cos V sinh v 


pCpS 


dt (sin V cos »+cosh v sinh v) 


Qz 


54=0 


Substituting these values for the constants of equation (Bl6) 
gives the equation for the velocity in dimensionless form: 


(B17) 


/ . vx 
/ sin — 

1 vx 
cosh — 

. 1 vx vx \ 

sink — cos — 1 

1 

r 

r r 1 

1 sin V 

1 cosh V 

sinh V cos v 1 

\sinh V 

COS V 

sin V ' cosh v / 


Differentiating equation (Bl7) and combining with 
equation (B7) result in an equation for the temperature in 
dimensionless form; 


tk 


tfk 

q"s 


i/ 

2v\ 


. VX T vx 
sin — cosh — 
r r 


. , vx vx 
sinh — cos — 
r r 


sinh V cos v sin v .cosh v 
sin V ~cosh v sinh v' cos v 


) 


(B18) 


The actual temperature cannot be determined from this 
equation, because the reference temperature t; is also a func- 
tion of 2. 

The difference between the wall temperature and the 
temperature at any point in the fluid at the same value of 2 
can be determined from equation (B 1 8) . From the definition 
of d*, equation (BI8) becomes 

e*= — 


ij 

i 457* 

f sinh V cos u-f-sinh— cos — 
r r , 

VX 

cosh V sin w-j- cosh — sin — I 
r r \ 

2v\ 

cosh V 

sin V ' 

sinh V 1 cos » | 


\ cos V 

‘sinh V 

sin V cosh v / 


(B19) 

The temperature difference wall-to-wall can be found from 
equation (Bl9) when x=r: 


! 

' kv \sinh v 


sin^ r-|-sinh^ v 


cosh r-f-sin v cos 




(B20) 


The temperature difference across the fluid for pure con- 
duction is 




(B21) 


The ratio of the actual temperature drop to the tempera- 
ture drop for pure conduction is 




dr 


{6. 


^ 

V)c V \sinh V 


sin^ y-f-sinh^ v 


cosh w-|-sin v cos v/ 


(B22) 


From the definitions of v and a. 
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